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Abstract 

Quantum computers are appealing for their ability to solve some tasks much faster than their 
classical counterparts. It was shown in [Aspuru-Guzik et al., Science 309, 1704 (2005)] that they, if 
available, would be able to perform the full configuration interaction (FCI) energy calculations with 
a polynomial scaling. This is in contrast to conventional computers where FCI scales exponentially. 
We have developed a code for simulation of quantum computers and implemented our version of 
the quantum full configuration interaction algorithm. We provide a detailed description of this 
algorithm and the results of the assessment of its performance on the four lowest lying electronic 
states of CH2 molecule. This molecule was chosen as a benchmark, since its two lowest lying 1 A\ 
states exhibit a multireference character at the equilibrium geometry. It has been shown that 
with a suitably chosen initial state of the quantum register, one is able to achieve the probability 
amplification regime of the iterative phase estimation algorithm even in this case. 
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I. INTRODUCTION 



Quantum chemical computations on conventional computers achieved a huge success in 
the last decades, both due to the enormous progress in computer technology and in the de- 
velopment of approximate quantum chemical methods. They became an indispensable part 
of basic as well as applied chemical research. The situation is however far from ideal. Exact 
solution of the Schrodinger equation within a finite one-particle basis set (full configuration 
interaction, FCI) is limited only to the smallest systems (diatomics, triatomics) because of 
the exponential scaling of its computational cost. Approximate methods with polynomial 
scaling must be used instead, of course at the cost of accuracy and reliability. 

Quantum computers on the other hand, if available, would offer exponential speedup for 
certain types of problems [JJ. The most prominent example is the famous Shor's algorithm 
[21 [3] for factoring integers, with potentially far-reaching consequences for cryptography. 
Another promising application of quantum computers, firstly suggested by Feynman [I], is 
efficient simulation of another quantum systems Abrams's and Lloyd's phase esti- 

mation method JTU] represents the polynomially scaling quantum algorithm for obtaining 
eigenvalues of local Hamiltonians. 

The first work connecting quantum computation and quantum chemistry by Lidar and 
Wang [11] concerned the efficient calculations of thermal rate constants of chemical reac- 
tions. This work in fact founded the new field of quantum chemistry: "Quantum chemistry 
on quantum computers" . Six years later, Aspuru-Guzik et al. [J_2] applied the Abrams's and 
Lloyd's algorithm [10] to quantum chemical energy calculations and proposed that quantum 
computers with 30 to 100 (noise free) qubits could already exceed the limits of classical FCI 
calculations. Since these two pioneering works, other papers involving energy calculations 
of excited states |13j . quantum chemical dynamics [H], calculations of molecular properties 
[To] , state preparation [TB],[T7] or global minima search [JS] were published. Whitfield et al. 
[T5] provided the detailed description of the quantum FCI algorithm [12] with emphasis on 
the factorization of an exponential of a Hamiltonian operator to elementary one and two 
qubit gates. Analysis of the resource requirements for the ground state energy calculations 
of the one dimensional transverse Ising model which took into account also the error correc- 
tion was done by Clark et al. [20] . An up to date summary of quantum chemical algorithms 
for quantum computers could be found in a recent review [2"T] . 
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Very recently, also the first proof-of-principle experiments representing quantum chemical 
calculations on quantum computers have appeared [221 123] ■ These are energy calculations 
of hydrogen molecule in a minimal basis set using optical quantum computer [22] arid NMR 
quantum computer |23j . In spite of the fact that these experiments treated only the simplest 
case of hydrogen molecule and did not have the resources to implement quantum error 
correction, it seems that simulations of physical phenomena on quantum computers might 
become available in foreseeable future |24|. 

The paper is organized as follows. The second section is devoted to the overview of 
quantum full configuration interaction algorithm, with the emphasis on the iterative phase 
estimation, since we feel that a more detailed description of the algorithm would be useful. 
The third section involves the application of our version of this algorithm to the four lowest 
lying electronic states of CH 2 molecule. 

II. OVERVIEW OF QUANTUM FULL CONFIGURATION INTERACTION 
(QFCI) ALGORITHM 

The first question to answer, when thinking about quantum chemical calculations on 
quantum computers, is how to map the quantum chemical wave function onto the register 
of quantum bits (qubits). Several different mappings have been proposed for this purpose 
[T2| IT3"] . Throughout this paper, we work with the simplest in terms of factorization of an 
exponential of the Hamiltonian to elementary gates and therefore the simplest for experi- 
mental realization, but the least economical one, so called direct mapping. In this approach, 
individual spin orbitals are directly assigned to qubits, since each spin orbital can be either 
occupied or unoccupied, corresponding to |1) or |0) states. Non-economical property lies in 
the fact that it actually maps the whole Fock space of the system (states with different num- 
ber of electrons) on the Hilbert space of the quantum computer. A compact mapping from a 
subspace of fixed-electron-number wave functions, spin-adapted or even symmetry-adapted 
wave functions to the qubits have been also proposed [121 [13], but the general factorization 
schemes for these mappings have not been discovered yet. 

Overview of the QFCI algorithm is divided into three subsections. In the first, state 
preparation is briefly discussed, the second is devoted to the phase estimation procedure 
and its improved iterative variant, which is the heart of the algorithm. The last section 
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shortly presents factorization of the exponential of the molecular electronic Hamiltonian to 
elementary gates. 

A. Initial state preparation 

Quantum FCI algorithm must be started with some initial guess state. Generally it holds 
that the closer is the initial guess to the exact wave function corresponding to calculated 
energy, the higher is the success probability of observing the energy. As is confirmed further 
in the paper, the Hartree-Fock guess, which is the first thing one tries to use, may not be 
successful especially in situations, where correlation (particularly the static one) plays an 
important role. In these situations, initial guess states from more sophisticated polynomially 
scaling methods can be used (e.g. CASSCF in a limited orbital space). 

Preparing a general initial state (vector from the Hilbert space of n qubits) is in fact a 
difficult task as this vector can contain up to 2 n non-zero components. Fortunately, initial 
guesses which include only a few determinants in a superposition (about 10 in our most 
difficult case) are sufficient for most purposes of quantum chemistry. These states can be 
prepared e.g. with the procedure described by Ortiz et al. [7] which scales as 0(N 2 ) 
in the number of determinants N . Preparation of general molecular-like states from the 
combinatorial space of dimension ( ) corresponding to distributing m electrons among n 
spin orbitals was presented in [17J. Preparation of many-particle states in a superposition 
on a lattice which can be then propagated by quantum chemical dynamics algorithm [14] 
was studied in [16]. 

Completely different approach is based on the adiabatic quantum computation [121 125] , 
which is an alternative to the classical gate model. Here one starts with a Hamiltonian, 
whose ground state is easy to construct, and gradually varies this Hamiltonian into the 
final exact one, where its ground state encodes the solution of the computational problem. 
If the variation of the Hamiltonian is slow enough, the system will remain in the ground 
state according to the adiabatic theorem. It has been shown that any quantum circuit 
can be simulated by an adiabatic quantum computation with an appropriately constructed 
Hamiltonian [2"B"l [2Tj . When preparing initial guess states for the quantum FCI algorithm, 
one starts with the Hartree-Fock Hamiltonian and evolves to the exact one [12]. This 
approach seems to be very promising especially due to the possible experimentally more 
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accessible realization. However, its scaling is rather unclear, because the maximum allowed 
speed of the adiabatic evolution depends on the minimum of the energy difference between 
the ground and the first excited state along the adiabatic path. 

B. Phase estimation 

The phase estimation algorithm is an efficient quantum algorithm for obtaining the eigen- 
value of an unitary operator U, based on a given initial guess of the corresponding eigen- 
vector. It is actually also the key part of the famous Shor's factoring algorithm (2J [3] . For a 
comprehensive description of the phase estimation procedure and citations to original papers 
see Nielsen and Chuang pp. 

Suppose that \u) is an eigenvector of U and that it holds 

U\u) = e 2 ^\u), 0G(O,1), (1) 

where <fi is the phase which is estimated by the algorithm. Quantum register is divided into 
two parts. The first part is the read-out part composed of m qubits on which the binary 
representation of (f) is measured at the end and which is initialized to the state |0)® m . The 
second part contains at the beginning of the procedure the corresponding eigenvector \u). 
The initial state of the quantum register thus reads 

Ireg) = |0) ® 10) ® ... ® |0)®|u), shortly |0)|u) (2) 

V v ' 

m qubits 

Application of Hadamard gates (tt/2 rotations) on all qubits in the first part of the register 
gives 

2 m -l 

l re s) = 7^E ( 3 ) 

After application of sequence of controlled-f/ 2 '" 1 operations (fcs are nonnegative integers 
from 1 to m), the register is transformed into 

keg) = ^=£ = ^=Y1 ^ ij4, \m- (4) 
vz 3=0 vz 7=0 
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The heart of the phase estimation algorithm is the inverse quantum Fourier transform 
(iQFT) [T] performed on the read-out part of the register which is transformed to |2 m 0)|ii). 
If the phase could be expressed exactly in m bits 



= 0.0102 • • • 0m (5) 
01 . 02 . . 0m , _ fn 1 , (a , 

it (and consequently the eigenvalue) could be recovered with unity probability by a mea- 
surement on the first part of the quantum register. 

If the desired eigenvector is not known explicitly (as is typically the case in quantum 
chemistry), we can start the algorithm with an arbitrary vector which can be expanded 
in terms of eigenvectors of U 



J2ci\ui). (7) 



The probability of obtaining 0j is due to linearity of the algorithm |cj| 2 . It is important to 
note that the initial guess does not influence the accuracy of the phase, only the probability 
with which the phase of a particular eigenstate is measured. 

The situation is more complicated when cannot be expressed exactly in m bits. Then 

= + 52~"\ (8) 

where = 0102 ■ ■ • m denotes the first m bits of the binary expansion and < 5 < 1 is a 
remainder. It can be shown (e.g. [28J) that the sum of success probabilities corresponding 

tO rounded down (P d own 0measured = 0) and rounded Up (P up 0mcasured = + 2 _m ) 

decreases monotonically for increasing m and in the limit m — >■ oo the lower bound is 

Pdown(5 = 1/2) + P up (5 = 1/2) = - + - > 0.81. (9) 

Note that linearity of the aforementioned scheme leads to the lower bound for -Pdown + -P U p 
corresponding to 0j being 0.81 • |q| 2 . 



6 



Iterative phase estimation (IPEA) 

With the use of semiclassical QFT j2H], the circuit representing the phase estimation 
algorithm can be greatly simplified, having only one ancillary qubit in the read-out part 
of the quantum register. The algorithm then proceeds in an iterative manner. The fc-th 
iteration of this scheme is presented in Figure [TJ 

The algorithm is iterated backwards from the least significant bits of <fi to the most 
significant ones, for k going from mto 1. The iteration again starts with the Hadamard gate 
on the read-out qubit followed by controlled-?/ 2 '' 1 operation. The equivalent of QFT is a 
single qubit z-rotation R z , whose angle Uk depends on the results of the previously measured 
bits 

RM) = (J e iL k ) (io) 

m— fe+1 , 

■* = - £ ^ (id 

followed by a Hadamard gate. 

IPEA is in fact completely equivalent to the original (multiqubit) phase estimation [28] . 
It thus suffers from the same decreasing of success probability when the phase cannot be 
expressed exactly in a particular number of bits. One possibility of a success probability 
amplification is performing more iteration steps (more than is the desired accuracy of 0): 
when extracting m' = m + log (2 + l/2e) bits, the phase is accurate to m binary digits with 
probability at least 1 — e pQ. This method is however not very useful, since implementing 
the U gate for large k is the algorithm's bottleneck in a realistic noisy environment [30J . 

Dobsicek et al. |30j came up with a different approach. They repeated the measurement 
for the least important bits of the phase binary expansion. Using the majority voting (for 
bit value or 1), the effective error probability decreases exponentially with the number of 
repetitions according to the binomial distribution. This measurement repetition only for the 
few least important bits of is unfortunately possible only if the exact eigenstates of U are 
available. 

When working with general initial states (the case of quantum chemistry), two scenarios 
presented in Figure [2] come into question. Maintaining the second part of the quantum 
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register during all iterations and amplification of the success probabilities by repeating the 
whole process is the first possibility. We denote this version as A version of IPEA. The 
biggest advantage of this approach is that one always ends up with one of the eigenstates of U 
in the second part of the quantum register as was the case in the original PEA. This happens 
through successive collapses of the system state into the corresponding eigensubspace. We 
demonstrate this approach in Figure [3] on the hydrogen molecule in a minimal basis set 
and random initial states (to avoid collapsing into the eigenstate which does not fit into 
the correct number of electrons - 2, we generated random initial states mixing only the 
components corresponding to correct number of electrons). The biggest disadvantage that 
might complicate its physical realization is the requirement for a long coherence time of the 
quantum register. 

Another possibility is to initialize the second part of the quantum register at every iter- 
ation step (B version of IPEA). Every iteration step (not only the least important bits of 
4> as in Ref. [3U]) must be repeated and measurement statistics performed. One could oth- 
erwise possibly mix bits belonging to different eigenvalues in different iterations and obtain 
an unphysical result. The biggest advantage of this approach is avoidance of the long co- 
herence times and therefore potentially easier physical implementation. On the other hand, 
the biggest disadvantage is that no improving of the overlap between the actual state of the 
quantum register and the exact wave function occurs during the iterations and one must 
"fight" the overlap problem at every iteration step. But as our simulations have shown and 
is discussed further, the situation concerning this overlap problem is in fact quite acceptable 
in practise. Small number of repetitions of each iteration is sufficient for amplification of the 
success probability to unity, when a suitable initial state of the quantum register is used. 

Application to molecular Hamiltonians 

The relevance of phase estimation for quantum simulations was first noticed by Abrams 
and Lloyd |10j . They took U in the form 

U = e irA , (12) 

where H is the Born-Oppenheimer electronic Hamiltonian and r is a suitable parameter 
which assures being in the interval (0, 1). 
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The electronic Hamiltonian can be expressed in the second quantized form as [31] 

1 , L 

pq pqrs X=l 

where h pq and (pq\rs) are one and two-electron integrals in the molecular spin orbital basis 
and a] and hi are fermionic creation and annihilation operators. The whole summation is 
formally expressed as a sum of individual terms hx- Since these creation and annihilation 
operators in general do not commute, exponential of the Hamiltonian cannot be written as 
a product of exponentials of individual hx, but a numerical approximation must be used 
[5]. The first-order Trotter approximation [32J can be expressed as 

L N 

e irH = e *rj: L x=1 ** = (J] e *xT/N\ + o(t 2 /N). (14) 

X=l 

By choosing N > (r 2 /e), we can implement U within an error tolerance of 0(e) using 
0(L(r 2 /e)) particular terms e thxT / N . Because L scales as 0(n 4 ) in the total number of spin 
orbitals n and individual terms e lhxT ' N can be built efficiently [7], the whole algorithm is 
efficient (polynomial). 

When taking into account that about twenty binary digits of <fi are necessary for quantum 
chemical applications, the maximal power of U in IPEA is 2 {rn -^ = 2 19 w 5 ■ 10 5 . This brines 
quite a big prefactor to the polynomial scaling of the algorithm and it is surely not necessary 
to emphasize that the Trotter approximation is the bottleneck of the algorithm. 

In our implementation, two external inputs are necessary. These are maximum (-E max ) 
and minimum (-E m i n ) energies expected in the studied system. We use U in the form 

U = e lT ( E —-"\ (15) 

where r reads 

(16) 



and the final energy is obtained according to the formula 



E = E max — . (17) 

T 
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This approach assures to be in the interval (0, 1). Aspuru-Guzik variant [12], where 
(f) = 0.5 for Hartree-Fock energy (E SC f), corresponds to the choice of E max = and E min = 
2-Escf- 

E min and E mSuX can be chosen arbitrarily, but one must be sure that the calculated energy 
is within this interval, otherwise one would end up with a non-physical energy, due to the 
periodicity of e 2m ^. The smaller the interval between the minimum and the maximum 
energies is, the less iterations of IPEA are necessary for desired precision of E and therefore 
less quantum gates are used. However, as the interval is smaller, r is bigger and more 



repetitions of the Trotter approximation (14) must be performed. In fact, when obtaining 
the energy with a fixed precision, the total exponential factor r • 2 m_1 is for the most 
"expensive" powers of U constant and independent of the size of the energy interval. 



Taking U in the form (15) does not pose any difficulties and indeed as the following circuit 



equality shows, just one more one-qubit rotation is needed. 



1 


p 




( 


i 




e ir(_B max -_f/) 




e -irH 








c 


girE max J 



C. Factorization of the exponential of a second quantized Hamiltonian 

Since we did not perform the factorization of the exponential of the molecular electronic 
Hamiltonian to one and two-qubit elementary gates in the numerical simulations, we will 
only briefly mention this topic, that was first studied by Ortiz et al. [7j. 

The first step is the Jordan- Wigner transformation [33] which maps the fermionic cre- 
ation and annihilation operators to spin operators represented by Pauli cr-matrices. The 
transformation has the form 



n-l 

n 



'n-l 



Gi cr 



+ ■ 



al \a n _. 



(18) 



n 

. 3=1 / \ 3=1 

where a± = l/2(a x ± ia y ) and the superscript denotes the qubit on which the matrix acts. 
After application of this transformation, the molecular Hamiltonian can be written in terms 
of strings of cr-matrices. It was shown that individual exponentials containing these strings 
can be built efficiently from one and two-qubit elementary gates [7J [191 E21 [31]. The circuit 
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representations of individual terms can be found in [19] or in the Supplementary Information 
of [22]- 

Once the factorization scheme is known, one can study the complexity of the whole 
algorithm. Scaling of the algorithm is given by the scaling of a single controlled action 
of the unitary operator ^ Ilx=i e t * lxT / N ^ from equation (14). Repetitions in the Trotter 
approximation (14) increase only the prefactor to the polynomial scaling, not the scaling 
itself. Also the required precision is limited, about twenty binary digits of 4> are sufficient for 
chemical accuracy. Detailed complexity analysis presented in [33] revealed that the algorithm 
scales as 0{n^, orh ) in the number of spin orbitals (number of qubits in the second part of the 
quantum register). This is in contrast to the exponential scaling of the full configuration 
interaction method on a conventional computer. 

The exponential speedup is demonstrated in Figure |4j Quantum gates count for a single 
controlled action of the unitary operator was done according to [19] and the number of gates 
corresponds to standard one and two-qubit gates from the set: Hadamard gate (one-qubit), 
CNOT gate (two-qubit), single qubit x-rotation R x (—ir/2), single qubit ^-rotation R z {6) 
and controlled R z {@) (two-qubit). We have to emphasize at this point that a simulation of 
the QFCI method on a conventional computer cannot achieve this exponential speedup in 
principle. Even though only polynomial number of elementary quantum gates is necessary, 
the simulation of an action of such single elementary gate scales exponentially itself, since 
it corresponds to a matrix vector multiplication of the dimension 2 n . 

Here we would like to note three more things. Firstly, we assumed that initial state 
preparation is an efficient step, as was already mentioned. Secondly, when a quantum 
chemical method with a scaling worse than 0(n 5 ) is used for calculation of an initial guess 
state on a conventional computer, then this classical step is becoming a rate determining 
one. Besides this, the classical computation of the integrals in the molecular orbital basis 
scales as 0(n 5 ) (due to the integral transformation) as well. 



III. APPLICATION TO FOUR LOWEST-LYING ELECTRONIC STATES OF CH 2 

Methylene molecule (CH 2 ) in a minimal basis set (STO-3G) is a simple, yet computa- 
tionally interesting system suitable for simulations and testing of the aforementioned QFCI 
algorithm. CH2 molecule is well known for the multireference character of its lowest-lying 
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singlet electronic state (a l A\) and is often used as a benchmark system for testing of newly 
developed computational methods (see e.g. [35 38J). When using the STO-3G basis set, 
the total number of molecular (spin)orbitals is 7(14). We therefore work with 15 qubits in 
the direct mapping approach (one qubit is needed in the read-out part of the register) . In 
spite of the fact that the complexity of simulations of the QFCI on conventional computers 
scales exponentially as the complexity of the classical FCI but with an order of magnitude 
larger prefactor [12], this system is still computationally feasible and for its properties an 
excellent candidate for one of the first benchmark simulations. 

Our aim was to verify the applicability of the QFCI for the ground as well as excited states 
exhibiting multireference character. We accordingly simulated the QFCI energy calculations 
of the four lowest-lying electronic states of CH2: X 3 B±, a l A%, b 1 B\, and c l A\. For CH2 
at the equilibrium geometry, the ground electronic state is not a closed-shell singlet, but a 
triplet state (X 3 Bi) with the electronic configuration 

(la 1 ) 2 (2a 1 ) 2 (16 2 ) 2 (3a 1 )(16 1 ). (19) 

The closed-shell singlet state (a l Ai), which can be qualitatively described by the electronic 
configuration 

(la 1 ) 2 (2a 1 ) 2 (l& 2 ) 2 (3a 1 ) 2 , (20) 

is the first excited state. This state exhibits the multireference character with the second 
important configuration 



(la 1 ) 2 (2a 1 ) 2 (16 2 ) 2 (16 1 



(21) 



The contribution of both closed-shell configurations becomes equal at linear geometries. The 



third electronic state (b x Bi) has the same spatial orbital configuration (19) as the ground 
state, but with singlet-coupled open shell electrons. The fourth electronic state (c 1 A 1 ) 
is represented by the same two configurations as the a 1 A± state but the amplitudes have 



the same sign and the amplitude of (21) is greater than that of (|l9|) [this state can be 



qualitatively described by the configuration (21)]. 

We simulated the QFCI energy calculations for C-H bond stretching (both C-H bonds 
were stretched, Figure pi), and H-C-H angle bending for a 1 Ai state (Figure pi). These 
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processes were chosen designedly because description of bond breaking is a hard task for 
many of computational methods and H-C-H angle bending since the a 1 A 1 state exhibits 
very strong multireference character at linear geometries. The equilibrium geometry of CH 2 
molecule was adopted from [39] and corresponded to r e = 1.1089 A and a e = 101.89°). 
Our work follows up the work by Wang et al. [13]. In this paper, the authors studied the 
influence of initial guesses on the performance of the quantum FCI method on two singlet 
states of water molecule across the bond-dissociation regime. They found out that the 
Hartree-Fock initial guess is not sufficient for bond dissociation and suggested the use of 
MCSCF method (CASSCF in particular). Few configuration state functions added to the 
initial guess improved the success probability dramatically. 

We also used and tested different initial guesses for QFCI calculations. Those denoted as 
HF guess were composed only from spin-adapted configurations which qualitatively describe 



certain state: in case of a A% configuration (20), in case of c A\ configuration (21), in 



case of X 3 Bi two triplet-coupled configurations (19) (with weights 1/2) and for b x Bi 
the same two configurations but singlet-coupled. Initial guesses denoted as CAS (a;,?/) guess 
were based on CASCI calculations with small complete active spaces (more details about the 
definition of the active spaces will be given further), which contained x electrons in y orbitals. 
Initial guesses based on CASSCF calculations as in [13] could be used in the same way. To 
be consistent, we employed the FCI wave functions in a limited active space composed 
of RHF orbitals, which were also used for the exponential of a Hamiltonian in the QFCI 
algorithm. Initial guesses were constructed only from the configurations whose absolute 
values of amplitudes were higher than 0.1. Those constructed from the configurations whose 
absolute values of amplitudes were higher than 0.2 are denoted as CAS(x,y), tresh. 0.2 
guess. All the initial guesses were normalized before the simulations. 

Similarly as in Ref. [12], the exponential of a Hamiltonian operator was implemented as 
a n-qubit gate. Factorization to the elementary one and two-qubit gates was performed only 
to examine the gate count, but not in the numerical simulations. We also did not take into 
account any decoherence and thus assumed that the exponential of the Hamiltonian can be 



obtained with an arbitrary precision by a proper number of repetitions in (14). One and 



two-electron integrals in the MO basis, parametrizing the Hamiltonian (13), were obtained 
using the restricted Hartree-Fock (RHF) orbitals. All ab initio calculations (FCI, RHF) 
were employed with our suite of quantum chemical programs 
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The phase in IPEA was always computed up to m = 20 binary digits. Maximum and 
minimum expected energies needed for the algorithm were set to E max = —37.5 a.u. and 
-Emm — —39.0 a.u. All presented success probabilities correspond to sum of the probabilities 
of rounding the phase up and down (P tot = P up + Pdown), therefore to probabilities of 
obtaining the final energy with precision « 1.43 • 10 -6 a.u. 

Finally, both of the aforementioned variants of IPEA (A and B) were tested. 

IV. RESULTS 

A. C-H bond stretching 

Results for the C-H bond stretching are summarized in Figures [7] - [9j 
Figure [7] presents the performance of the A version of IPEA with maintaining the second 
part of the quantum register during all iterations. Subfigures a - d represent the simulations 
of the energy calculations of the four electronic states: a: a l A\, b: c l A\, c: X 3 Bi, and d: 
b 1 B\. Overlap between the initial HF guess wave function and the exact FCI wave function 
as well as this overlap scaled by the factor 0.81 [according to ([9])] are shown. Figure [7] also 
presents the success probabilities of IPEA for the HF initial guess and initial guesses based 
on the CASCI calculations with certain small complete active spaces. Definition of these 
active spaces is complicated by the fact that swapping of molecular orbitals occurs when the 
C-H bonds are prolonged. To maximize the overlap between the initial and the exact wave 
functions, we constructed the active spaces from the actual highest occupied and the lowest 
unoccupied molecular orbitals at a given geometry. For X 3 P 1 , b 1 Bi and c 1 v4 1 , where 
the l&i orbital is involved in the qualitative description of the state (at the equilibrium 
geometry), this orbital was always included in the active space [l&i orbital which is the 
LUMO (5th molecular orbital) at the equilibrium geometry becomes the 7th when going to 
three times prolonged C-H bonds]. Definition of the complete active spaces is summarized 
in Table |TJ Dotted line corresponding to the probability 0.5 bounds the region where the 
algorithm can be safely used and the total probability amplified by repeating the whole 
process. 

Figures [8] and [9] present the performance of the B version of IPEA. In this version, the 
second part of the quantum register is reinitialized at every iteration step. Figures [8] and 
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state 




p A C//| A \ 
LAS^4,4 ) 




a l Ax 


highest occupied MO, 


two highest occupied MOs, 




lowest unoccupied MO 


two lowest unoccupied MOs 






highest occupied MO, 


two highest occupied MOs, 




c 1 A 1 


unoccupied l&i MO 


unoccupied l&i MO, 






lowest unoccupied MO 
(other than lb±) 








two highest occupied MOs, 


two highest occupied MOs, 


X 3 5i 




unoccupied lb\ MO, 


three lowest unoccupied MOs 




lowest unoccupied MO 
(other than lb±) 


(including l&i) 






two highest occupied MOs, 


two highest occupied MOs, 


b 1 B 1 




unoccupied 16i MO, 


three lowest unoccupied MOs 




lowest unoccupied MO 
(other than lb±) 


(including l&x) 



Table I: Summary of the complete active spaces (CAS) used for the calculations of initial guesses 
for IPEA (Figures [7] and [9]), occupation/unoccupation refers to the lowest closed-shell configuration 



(20). 



[9] demonstrate the success probabilities for different number of such repetitions (11-101). 
Figure [8] shows the results and limits of the HF guess for a 1 Ai state. Figure [9] presents the 
results of the "best" initial guesses in terms of price/performance ratio for all four states. 



B. H-C-H angle bending 



Results for the H-C-H angle bending are summarized in Figures 10 and 11 Simulations 
concerning this process involve only the a l A\ state as this state exhibits a strong multire- 
ference character when going to linear geometries. In this case, no swapping of molecular 
orbitals occurs during the process and the complete active space CAS(2,2) was always con- 
structed from 3ai (HOMO) and l&i (LUMO) molecular orbitals. Moreover, due to the 
different symmetry of these orbitals, only two configurations contribute to CAS(2,2) wave 



function: doubly occupied HOMO [configuration (20)] and doubly occupied LUMO [config- 



uration (21)]. Both of these configurations have for all values of a (H-C-H angle) absolute 
values of amplitudes higher than 0.2. 



Figure 10 presents the results of the A version of IPEA. Overlap and scaled overlap of 
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the initial HF guess wave function and the exact FCI wave function is again shown as well as 
the success probabilities for HF and CAS(2,2), tresh. 0.2 guesses and dotted line bounding 
the safe region. Performance of the B version of IPEA with HF and CAS(2,2), tresh. 0.2 



initial guesses is illustrated in Figure 11 



V. DISCUSSION 
IPEA - A version 



Results of the simulations with A version of IPEA numerically confirm that success 
probabilities always lie in the interval | (VWt Inexact) | 2 ■ (0.81, l), depending on the value of 
the remainder 5 This algorithm can be safely used when the resulting success probability 
is higher than 0.5 (as it can then be amplified by repeating the whole process). We would 
like to note that when studying the applicability of this algorithm, one must monitor the 
scaled overlap between the initial guess and the exact wave function. 

Success probability higher than 0.5 is securely fulfilled with the HF initial guess at the 
equilibrium geometry for all four simulated states. When going to more stretched C-H bonds 
or linear geometry, RHF initial guess fails. CAS (2,2) initial guess improves the success 
probability in case of a l Ai and c 1 Ai states near the equilibrium geometry but in the region 
of more stretched C-H bonds it also fails. In this region, CAS (4,4) initial guesses must 
be used. For a l A\ state, CAS(4,4), tresh. 0.2 guess is sufficient but for c 1 A%, even the 
CAS(4,4) guess fails for few lengths of C-H bonds: r/ro = 2.2,2.8 — 3.0. For these points 
of the potential energy surface, where the overlap between the initial guess and the exact 
wave functions is not high enough, bigger active space should probably be used. 

The situation is more difficult for the states of Bi symmetry (X 3 Bi, b l Bi) when the 
C-H bonds are stretched. Here even CAS (4,4) initial guess fails and bigger active space - 
CAS (4,5) - must be used for initial guess state calculations. In STO-3G basis set, this bigger 
complete active space contains five from the total number of seven molecular orbitals and 
represents therefore nearly the whole space. For this reason, we performed the classical FCI 
calculations with the cc-pVDZ basis set (Is orbital on the carbon atom was kept frozen to 
reduce computational demands), where the total number of molecular orbitals is 24 (much 
more than the number of molecular orbitals in the complete active space), and verified 
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that the overlap between the CAS (4,5) and the exact wave function is sufficiently high, 
essentially the same as in STO-3G basis set. Apart from the active space size, initial guess 
states always contained at most 12 configurations, but usually 8 or even less for nearly 
dissociated molecule. This observation is in agreement with the results of [13J, where few 
configuration state functions added to the initial guess improved the success probability 
dramatically. 

High success probabilities (over 0.8) can on the other hand be obtained with CAS(2,2), 
tresh. 0.2 initial guess for a 1 Ai state during H-C-H angle bending. Initial guess states for 
this process correspond to only two configurations and are thus very easy to prepare (e.g. 
according to [7j). 

We would like to note one more thing about the A version of IPEA. This version can 
be in principle effectively used for ground state energy calculations even when the success 
probability is lower than 0.5. In fact, since it always converges to one of the eigenstates, it 
can be used in all situations when the overlap between the initial guess state and the exact 
wave function is not exponentially small (with respect to size of the system). In those cases, 
one will repeat the procedure (a polynomial number of times, depending on the value of 
the aforementioned overlap) and look for the lowest energy. The only thing one must pay 
attention to is to set the limits E max and £ m ; n so that all the eigenvalues of the Hamiltonian 
lie between them. Then only it is sure that the lowest energy is not artificial and really 
corresponds to the lowest eigenvalue. 

IPEA - B version 

This version of IPEA is characteristic by repeated initial state preparation in each iter- 
ation and has the disadvantage that no collapsing of the system and improving the overlap 
between the actual state of the quantum register and the exact wave function occurs. The 
situation is however not so bad because one does not "fight" against the overlap between 
the initial guess and unwanted eigenfunctions at every iteration. This would happen only 
if all binary digits of the phase were opposite to binary digits of the phases of all other 
eigenstates, which is clearly not possible. Our simulations proved that a relatively small 
number of repetitions (~ 51) at each iteration step give nearly unity success probability 
when modest-size initial guesses are used. 
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,) nicely show their limits (for a l A\ 
state). These are: r/r = 2.3 for C-H bond stretching and 170° for H-C-H angle bending. 
We have chosen the "best" initial guess states in terms of price performance ratio for each 
of the four electronic states: a 1 A\ (C-H bond stretching), c 1 A\\ CAS(4,4), tresh. 0.2; 
X 3 Bi, b 1 B l : CAS(4,5), tresh 0.2; a l A 1 (H-C-H angle bending): CAS(2,2), tresh. 0.2. 
They were chosen to contain the minimum number of configurations, yet give high enough 
success probabilities. These initial guesses performed very well and the only exception are 
a few points for c 1 Ai state, where a bigger active space is desirable. 

We have not simulated the decoherence phenomena and when taking it into account the 
situation will surely change. Quantum error correction [4T] would probably be needed for 
the A version, which would increase the number of required qubits as well as the complexity 
of the quantum circuit, while B version should be more robust. 

VI. CONCLUSIONS 

We have developed programs (C++) for simulation of a quantum computer and success- 
fully performed QFCI energy calculations of ground as well as excited states of CH 2 molecule 
that exhibits multireference character. We have demonstrated that energies at equilibrium 
geometry are accessible with RHF initial guesses, which are easy to prepare. CASCI initial 
guess states with small complete active spaces composed of relatively few configurations 
(« 10) are sufficient even for a nearly dissociated molecule to achieve the probability ampli- 
fication regime of the IPEA algorithm. Its version B with repeated initial state preparation 
seems to be a better candidate for the first real larger-scale QFCI calculations than the A 
version, since it does not require a long coherence time. 



The results with the RHF guesses (Figures M and 
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Figure 1: The fc-th iteration of the iterative phase estimation algorithm (IPEA). The feedback 
angle depends on the previously measured bits. Note that k is iterated backwards from into 1. 
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a) Maintaining the second part of the quantum register during all iterations. 
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b) Initial state preparation per iteration. 



Figure 2: Comparison of the two versions of IPEA. Oracle | ISP | prepares the initial state {Initial 
State Preparation). 
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Figure 3: Energies of the four electronic states of H2 in STO-3G basis which were obtained by 
QFCI (IPEA version A) with randomly generated initial guess states. Small figure inside presents 
the increasing overlap between the actual state of the second part of the quantum register and the 
exact wave function for one of random runs of the algorithm leading to the ground state. 
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Figure 4: The exponential speedup of the QFCI over the FCI. In case of the FCI (blue), dependence 
of the number of Slater determinants in the FCI expansion on the number of basis functions is 
shown. In case of the QFCI (red), dependence of the number of one and two-qubit gates needed 
for a single controlled action of the unitary operator (for details see the text) on the number of 
basis functions is presented. There are logarithmic scales on y axes and the points in the graph 
correspond to depicted molecules (hydrogen, methylene, methane, ethane, and benzene) in the 
cc-pVDZ basis set. 
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Figure 5: Energies of the four simulated states of CH2 for the C-H bond stretching, tq denotes the 
equilibrium bond distance. 
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Figure 7: Success probabilities of the A version of IPEA for the four electronic states of CH2 and 
different initial guesses, tresh 0.2 means that only configurations with absolute values of amplitudes 
higher than 0.2 were involved in the initial guess, ro denotes the equilibrium bond distance. 
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Figure 8: Success probabilities of the B version of IPEA with HF guess for a A\ state and different 
number of repetitions of individual bit measurements, tq denotes the equilibrium bond distance. 
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Figure 9: Success probabilities of the B version of IPEA with "best" initial guesses and different 
number of repetitions of individual bit measurements for all four states, tq denotes the equilibrium 
bond distance. 
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Figure 10: Success probabilities of the A version of IPEA for the a 1 A\ state with HF and CAS(2,2) 
initial guesses, tresh 0.2 means that only configurations with absolute values of amplitudes higher 
than 0.2 were involved in the initial guess, a denotes the H-C-H angle. 
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Figure 11: Success probabilities of the B version of IPEA with HF and CAS(2,2), tresh. 0.2 initial 
guesses and different number of repetitions of individual bit measurements for the a l A\ state, a 
denotes the H-C-H angle. 
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